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A survey of pairing properties of nucleonic matter is presented that includes the off-shell propa- 
gation associated with short-range and tensor correlations. For this purpose, the gap equation has 
been solved in its most general form employing the complete energy and momentum dependence of 
the normal self-energy contributions. The latter correlations include the self-consistent calculation 
of the nucleon self-energy that is generated by the summation of ladder diagrams. This treatment 
preserves the conservation of particle number unlike approaches in which the self-energy is based on 
l/^ , the Brueckner-Hartree-Fock approximation. A huge reduction in the strength as well as temperature 

■ and density range of ^Si-'^Di pairing is obtained for nuclear matter as compared to the standard 

' BCS treatment. Similar dramatic results pertain to pairing of neutrons in neutron matter. 
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I. INTRODUCTION 



Advances in the understanding of the single-particle (sp) properties of nucleons in nuclei and nuclear matter [J 
demonstrate the dominant influence of short-range and tensor correlations in generating the distribution of the spectral 
• strength. One aspect of this influence is expressed in the global depletion of Fermi sea due to these correlations. A 
' recent experiment from NIKHEF puts this depletion of the proton Fermi sea in ^°^Pb at a little less than 20% |3| in 
CO 1 accordance with earlier nuclear matter calculations 0] . Another consequence of the presence of short-range and tensor 
correlations is the appearance of high-momentum components in the ground state to compensate for the depleted 
, strength of the mean field. Recent JLab experiments Q indicate that the amount and location of this strength is 

■ consistent with earlier predictions for finite nuclei |^ and calculations for infinite matter ■ Of particular relevance 
' is the observed energy distribution of the sp strength of the high-momentum nucleons that is located at energies far 

■~>«^ ' removed from the Fermi energy. This situation is similar to the case of nuclear matter |^ and also holds above 
1-^, ' the Fermi energy. Such sp strength distributions lead to substantial modifications of the calculated properties of 
I , . low-energy phenomena, like pairing, as compared to those generated by a traditional mean-field treatment. 
^ ' Pairing properties of nuclear and neutron matter have been studied for quite some time. A recent review can be 
^ , found in Ref. |9j . The study of pairing correlations in infinite quantum systems is an interesting question in general 

■ relating nuclear and neutron matter to the classical systems studied in condensed matter physics as well as to the 
^ ' studies of pairing correlations in fermionic atomic gases |lOl | . The study of pairing correlations in nuclear and neutron 

• '~j matter is of particular interest for the understanding of properties of neutron stars. The formation of a BCS gap has 
rS a significant effect on the neutrino emissivity, which is crucial for the cooling of neutron stars pdj . The existence of 

. superfluid layers will affect the rotation of neutron stars. 
" " ' Of particular interest is the observation that several calculations for the bare nucleon-nucleon (NN) interaction in 
the ^Si-^Di coupled channel lead to a sizable gap of around 10 MeV at normal nuclear matter density 0, 0, 0, 
Those calculations typically solve the BCS gap equation using a mean- field sp propagator with sp energies 
as determined e.g. in a Brueckner-Hartree-Fock (BHF) calculation. 

The empirical data of finite nuclei do not exhibit any indications for such strong proton-neutron pairing correlations, 
which would correspond to a gap as large as 10 MeV. Therefore one may conclude that such evaluations of a pairing 
gap in infinite matter yield quite different results from those observed for finite nuclei. This appears plausible since 
pairing effects are very sensitive to the sp spectrum close to the Fermi energy, which is continuous in infinite matter 
while it is discrete in finite nuclei. 

This observation may also cast some doubt on the reliability of the approach to determine the pairing gap for infinite 
matter as outlined above. Therefore attempts have been made to go beyond this mean-field approach. One issue 
discussed in the literature has been the role of vertex corrections, i.e. the medium dependence of the NN interaction 
to be used in solving the gap equation 0,0,0,1131. Using effective interactions like the Gogny force 21] Shen 
et al. 0, ll^l observe indeed a significant effect, which is larger for nuclear matter than for neutron matter. The 
effects of the so-called induced interaction could indeed affect the low-energy spectroscopy of nuclear matter in a 
significant way. Before conclusions can be drawn, calculations employing realistic interactions should be performed, 
which, unfortunately, are very difficult |22j |. In addition, it is possible that polarization effects are different in nuclear 
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matter and finite nuclei on account of the difference of the sp spectra discussed above. 

In the present work we would like to focus the attention to the study of pairing correlations with a proper treatment of 
the nucleon propagator that accounts for the effects of short-range and tensor correlations. The crucial effects of short- 
range and tensor correlations on the distribution of the sp strength have recently been treated at a very sophisticated 
level. Various groups have developed techniques to evaluate the sp strength from realistic NN interactions, like the 
Argonne V18 or the CDBonn interaction, within the self-consistent Green's function (SCGF) approach 
I25I l26l I27I I2M I30I I31I Is^ l . Such calculations reproduce the energy- and momentum-distribution of the sp strength 
corresponding to high-momentum nucleons as observed in experiments j3| and account for the depletion of the mean- 
field strength obtained in . 

In these calculations the scattering equations for two nucleons in the medium are solved employing dressed sp 
propagators that contain the complete information about the energy- and momentum-distribution of the sp strength. 
The resulting scattering matrix T is then used to calculate the nucleon self-energy. The solution of the Dyson equation, 
employing this complex and energy dependent self-energy, yields the sp propagator that enters the T matrix equation 
to close the self-consistency cycle. In determining the two-nucleon Green's function or the corresponding scattering 
matrix T, one has to face the problem of the so-called pairing instabilities that reflect the existence of NN bound 
state solutions in the scattering equation. The presence of such solutions provides a major numerical obstacle for a 
self-consistent evaluation of one- and two-body propagator within the normal Green's function approach at densities 
where such instabilities can occur. For this reason, recent SCGF calculations have been performed for temperatures 
above the critical temperature for a possible phase transition to pairing condensation ^ .25j . 

In Sec. II of this work we will review some basic features of the SCGF method for the normal Green's function at 
temperatures above the critical one for the pairing instability. The effects of the anomalous Green's functions will 
be considered in Sec^ III. This leads to the generalized gap equation that corresponds to the homogeneous scattering 
equation for dressed nucleons. The analogous problem occurs for two nucleons in the vacuum, where the solution 
of the homogeneous scattering equation yields the description of the bound two-nucleon state, the deuteron. In the 
nuclear medium, however, one has to consider a scattering equation with dressed sp propagators. If one considers the 
HF or BHF approximation for the normal sp propagator, in which the spectral function is replaced by a (5-function, 
this gap equation reduces to the usual BCS approximation. 

Results of SCGF calculations above the critical temperature will be discussed in Sec. IV. We will pay special 
attention to the spectral functions associated with sp strength around the Fermi energy and discuss, to which extend 
these exhibit typical precursor phenomena for a phase transition to a pairing condensate [33.] . In Sect. IV we will 
also investigate how the distribution of the sp strength modifies the solution of the gap equation and the subsequent 
predictions for a phase transition to a superfluid state of nuclear matter. Such investigations have been performed 
before by Bozek 34, 35] using simplified models for the NN interaction. These studies suggest an important sensitivity 
to these correlations, accompanied by a substantial suppression of the strength of pairing. A similar conclusion was 
reached for a realistic interaction in Ref. [s^ by studying the phase shifts of dressed nucleons in the medium that 
signal the presence of bound pair states as in the case of two free particles. The inclusion of dressing effects in the 
study of pairing has also been studied in Refs. Is^,!!!,!!^ based on the hole-line expansion of the nucleon self-energy. 
While also in this work a substantial reduction of the strength of pairing is observed, the implementation of the 
scheme to solve the gap equation relies on approximations that do not conserve particle number, since they involve 
the introduction of quasiparticle strength factors to represent the effect of dressing. Final conclusions are drawn in 
Sec. V. 



II. GREEN'S FUNCTIONS AND T-MATRIX APPROXIMATION 



One of the key quantities within the SCGF approach is the sp Green's function, which can be defined in a grand- 
canonical formulation for both real and imaginary times t, t' |40l ]: 

iCU t- x' t') = tr{cxp[-/3(g - ^^N)m{^t)4>H^'m 

^ ' ' ' tT{exp[-f3{H - fiN)]} ' ^ ^ 

where T is the time ordering operator, (3 the inverse temperature, and fi the chemical potential of the system. 
Due to the invariance of the trace under cyclic permutations the one-particle Green's function obeys the following 
quasi-periodicity condition 

G(x, t = 0; x', <') = -e'^^G(x, t = -i/3; x', t'). (2) 



For a system invariant under translation in space and time, the propagator depends only on the differences ]x — x'] and 
t — t', allowing a Fourier transformation into momentum and energy variables k and lo. Due to the quasi-periodicity 
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in Eq. the Green's function can be expressed in terms of the Fourier coefficients G{k, z^), where = are the 
(fcrmion) Matsubara frequencies with odd integers ly. Since these are related to the spectral function A{k,uj) by 

G(k,z^) = / 77- 77, (3) 

G{k, Zu) can be continued analytically to all non-real z. On the other hand, the spectral function is related to the 
imaginary part of the retarded propagator G{k, uj + irf) by 

A(fc,w) = -2ImG(fc,w + i?7). (4) 

In the limit of the mean-field or quasi-particlc approximation the spectral function is represented by a J-function and 
takes the simple form 

A{k, uj) = 2ii5(uj + ^ - Efc) = 2tt5{uj - Xfe) , (5) 

with the quasi-particle energy for a particle with momentum k and Xfc = £fe ~ Note, that for convenience we 
define the energy variable relative to the chemical potential /i. 

The sp Green's function is obtained as a solution of the Dyson equation, which, for a translationally invariant 
system, is a simple algebraic equation 
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k 

UJ + ^ - E(fc,Lj) 

2m 



G(fc,c.) = l, (6) 



where w) denotes the complex self-energy. By expanding the self-energy in terms of one-particle Green's functions 
it can be demonstrated that it inherits all analytic properties of G. It it thus possible to write 

m = ^"^{k) - 1 du;' ^'"'^(^^' + '"^ . (7) 

The next step is to obtain the self energy in terms of the in-medium two-body scattering T matrix. It is possible to 
express ImE(fc, lu + iry) in terms of the retarded T matrix |^l4Cl|| 

ImE(fc,c. + ir,) ^ - J J — (kk'|Imr(c. + c.' + i,/)|kk') 

x[f{Lo') + b{uj + uj')]A{k\u;'). (8) 



Here and in the following 



e 



m = (9) 

denote the Fermi and Bose distribution functions, respectively. The pole in the Bose function 6(f2) at $1 = is 
compensated by a corresponding zero in the T matrix |4ll |42| such that the integrand remains finite as long as the T 
matrix does not acquire a pole at this energy. Such a pole may occur below a critical temperature Tc, a phenomenon 
that is often referred to as a pairing instability. We will come back to this problem below. 
The scattering matrix T is to be determined as a solution of the integral equation 

(kk'|r(f]-|-i,/)|pp') = (kk'll^lpp') + / ^!l^(kk'l^|qq')Gfi(qq',l] + *77) 

J (27rj 

x(qq'|T(17-Mr/)|pp') , (10) 

where 

J-00 27r 27r il - uj - u' + it] 



stands for the two-particle Green's function of two non-interacting but dressed nucleons. 
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The Green's function method yields a hierarchy of relations for the A''-particle Green's functions. The Dyson 
equation for the one-particle Green's function involves the two-body potential as well as the two-particle Green's 
function Gii(x, i; ■ • • ;x"',i"'). In general, the equation of motion for the A^-particle propagator will be coupled to 
the (N + l)-particle propagator, if the Hamiltonian contains a two-body interaction. In the self-consistent T-matrix 
approach, as outlined above, one ignores the effects of A^-particle Green's functions with N larger equal to three, but 
solves the coupled equations for the one- and two-body Green's functions in a self-consistent way. 

In order to allow for an efficient solution of the two-body scattering equation, one expresses the two-particle Green's 
function in Hll|l as a function of the total momentum P = ^(ki -I- and the relative momentum q = ^(ki — k2). 
Using the usual angle-average approximation for the angle between P and q (see e.g. ji^l for the accuracy of this 
approximation), the two-particle Green's function can be written as a function of the length of these two- vectors, P 
and q, only. This approximation leads to a decoupling of partial waves with different total angular momentum J. 
Therefore the integral equation IjlUI) reduces to an integral equation in only one dimension of the form 

9 r°° 

{q\T,i,'^{P,n + irj)\q') = {q\V,fT^q') + J2^ [ dk' k'^ {q\Vd,F \k') G^iP^n + i^^k') 

x(fc'|T/,f,^(P,17 + i,y)|g'). (12) 

The summation of the partial waves, 

(kk'|lmT(r! + i77)|kk')-^ ^ (2J+l)(2T+l)(g(k,k')|ImT,f^(P(k,k'),r! + i?7)|g(k,k')), (13) 

(JST)l 

yields the T matrix in the form that is needed in Eq. ©. Finally, the Hartree-Fock contribution has to be added to 
the real part of S 

^"^i^)-^ E (2J+1)(2T+1) / i^(g(k,k')|l^f^^k(k,k'))n(fc'), (14) 

(.]ST)l ^ ' 

where n{k) is the correlated momentum distribution 

/ + 00 1 
f^/(c.)A(fc,c.). (15) 

Note, that Eq. H14|l corresponds to a generalized Hartree-Fock contribution, since the full one-particle spectral function 
is employed. 



III. PAIRING IN THE T-MATRIX APPROXIMATION 



For temperatures below the critical temperature for a transition to a superfluid one has to supplement the evaluation 
of the normal Green's function G{k,ijj) with the anomalous Green's function F{k,uj). While the self-consistent 
inclusion of ladder diagrams has reached quite a sophistication, it remains to fully account for the possibility of 
a pairing solution in such calculations for realistic NN interactions. In this paper we present a first step towards 
such a complete scheme by including the full self-consistent dressing due to normal self-energy terms generated by 
ladder diagrams, in the calculation of the anomalous self-energy and the solution of the corresponding generalized 
gap equation. 

The inclusion of the anomalous Green's function F(k, lu) yields a modification of the normal Green's function in 
the superfluid phase that can be written as [4J, |4^ |43, |43 

Gs{Kuj + ir]) = G{k,u} + ir]) - G{k,u} + ir])A{k)F{k,uj + ir]) 

F{k,uj + ir]) = G{-k,~uj~ir])Gs{k,uj + irj)A{k). (16) 

under the assumption that the anomalous part of self-energy A does not depend on the energy. Therefore the full 
Green's function can be obtained as 

^'^^^ " + ^''^ = G(k,c. + ^,y)-l + A^2(k)G(-k,-c.-^,) ' ^^^^ 
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These equations must be supplemented with the definition of the anomalous self-energy 

^(p) ^ / £ / (0^ ^P'^'^^ ^^"^^^^^ + "'^^^'^^ ■ ^^^^ 

If one employs Eq. pt)|l and uses the representation of the Green's function G in terms of the spectral function in 
Eqs. Q and supplemented by a corresponding definition of a spectral function for the total Green's function 

As{k,uj)^-2lmGs{k,u} + iTj), (19) 

the expression for the self-energy A can be rewritten j3^] in a partial wave expansion 

^ ^ Jo J-oo 27r J^^ 2Tr ' -uj - uj' 

If we ignore for a moment the difference between the spectral functions A and A^, we see that this equation for the 
self-energy A corresponds to the homogeneous scattering equation for the T-niatrix in (|10|l at energy = and 
center-of-mass momentum P = 0. This means that a non-trivial solution of Eq. (|20|l is obtained if and only if the 
scattering matrix T generates a pole at energy 51 = 0, which reflects a bound two-particle state. This is precisely 
the condition for the pairing instability discussed above, demonstrating that this treatment of pairing correlations is 
compatible with the T-matrix approximation in the non-superfiuid regime discussed in Sec. II. 

We may also consider Eq. H20() in the limit in which we approximate the spectral functions A{k,u}) and As{k,uj) 
by the corresponding mean-field and BCS approximation. The expression for the normal spectral function has been 
presented already in Eq. |(SJ|. The BCS approximation for the spectral function yields 



Mk, ^) - 2^ (^^IJ^^C'^ - ^^0 + + ) . (21) 



with the quasi-particle energy 



E^ = Jxl + A^{k). (22) 



Inserting these approximations for the spectral function into Eq. (|20|l and taking the limit T = reduces to the usual 
BCS gap equation 



■■"^{p)=Y.- rdkk^p\v,r\k)^A-r{k). (23) 



Therefore we can consider Eq. (|20|l as a generalization of the usual gap equation. It accounts for the spreading of sp 
strength leading to a generalization of the form 



/ TT / -T^A{k,Lu)A,{k,Lu')- 



,....,„.,.....„y)^-^(-)-^(-r (24) 

-2Ek J-oo 27r J_^ 2tt -lo - uj' 



IV. RESULTS AND DISCUSSION 
A. SCGF above the critical temperature 



In the first part of this section we will focus the attention to the discussion of SCGF calculations for symmetric 
nuclear matter and pure neutron matter at temperatures above the critical temperature Tc for a phase transition to a 
pairing condensate. This means that we solve Eqs. (|12|l - (|15|l in an iterative procedure until a self-consistent solution 
is obtained. Some details of this procedure have been published in @ and |43 |. 

As a typical example we present in Fig. ^ results for the imaginary part of the retarded self-energy Ti{k,ui) for 
nucleons with a fixed momentum fc as a function of the energy variable uj. All the results displayed in this figure have 
been determined for symmetric nuclear matter at the empirical saturation density using the CDBonn |24j interaction. 
The energy scale in this figure has been constrained to energies around the Fermi energy (w = 0), since the self-energy 
is most sensitive to the temperature for these energies. The results for temperatures larger than or equal to 4 MeV, 
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FIG. 1: Imaginary part of the retarded self-energy for nucleons with momentum k = 225 MeV/c in symmetric nuclear matter 
at the empirical saturation density (p — 0.16 fm~'^). Results for temperatures T larger than or equal to 4 MeV are directly 
determined by SCGF calculations, while the T = result originates from an extrapolation. All results were obtained using the 
CDBonn interaction. 



which are all above Tc (see below), have been obtained directly from SCGF calculations. They exhibit a rather smooth 
dependence on the temperature, so that an extrapolation to temperatures below Tc appears feasible. As an example 
of such an extrapolation we show the T = result. This extrapolation has been done with the constraint that the 
imaginary part of the self-energy for T — vanishes at lo =0. 

The imaginary part of the nucleon self-energy is also displayed in the upper panel of Fig. |21 The purpose of this 
figure is to visualize some differences between various models of the NN interaction and between symmetric nuclear 
matter and pure neutron matter. Therefore we consider a larger interval for the energy variable uj. The imaginary 
parts of the self-energy derived from the CDBonn interaction and the Argonne V18 (ArV18) interaction 23] are very 
similar at energies around oj — 0. At those energies the ArV18 yields a slightly weaker imaginary part than CDBonn. 
The differences get larger at positive values of uj, where the imaginary part of the self-energy derived from ArVlB 
reaches a minimum of around -100 MeV at an energy lo around 1.7 GeV. The minimum for the CDBonn interaction 
is only about -35 MeV and occurs at energies uj around 0.5 GeV. 

This is another indication of the feature, that a fit of a local interaction, like ArVlS, to NN phase shifts yields a 
larger amount of NN correlations than a fit of a non-local relativistic meson exchange model like the CDBonn fitting 
identical phase shifts [iol f50l|. In short: The ArVlS is a stiffer interaction than CDBonn. A further illustration is 
provided by a Hartree-Fock calculation for nuclear matter at the empirical saturation density that yields a total energy 
30 MeV per nucleon for the ArVlB, while CDBonn generates 5 MeV per nucleon This implies that the generalized 
Hartree-Fock contribution to the self-energy in Eq. {Tj), defined in Eq. {71), is more repulsive for ArVlB, and a larger 
part of the attraction is provided by the energy-dependent contribution to the real part of the self-energy. This is 
immediately obvious, since the energy-dependent contribution to the real part of S is connected to the imaginary part 
by a dispersion relation. The lower panel of Fig. [21 which displays results for the real part of the self-energy, illustrates 
this observation: The energy dependence is larger for ArVlS as compared to CDBonn. The weaker attraction of the 
self-energy derived from CDBonn (for most values of uj) reflects the less repulsive contribution of the generalized 
Hartree-Fock contribution. 

Figure|21also displays results for the real and imaginary part of the self-energy for neutrons with the same momentum 
{k = 225 MeV) in pure neutron matter. The density of neutron matter considered in this figure is one half of the 
empirical saturation density of nuclear matter, which implies that these systems have the same Fermi momentum. 
The imaginary part of the self-energy in neutron matter is weaker than the corresponding one for symmetric nuclear 
matter, reflecting the dominance of proton-neutron correlations. For both interactions a minimum is obtained around 
1.7 GeV. At these high energies the absolute value for the imaginary part of the self-energy is about a factor three 
larger for ArVlS than for CDBonn. This means that the distribution of sp strength to high energies due to central 
short-range correlations is much stronger for the local ArV18 interaction than for the non-local meson-exchange model 
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FIG. 2: Imaginary (upper panel) and real part (lower panel) of the retarded self-energy for nucleons with momentum k = 
225 MeV/c in symmetric nuclear matter at the empirical saturation density (p = 0.16 fm~^). Results obtained for CDBonn 
interaction are compared to those resulting from SCGF calculations using ArV18. Also included are results for neutrons with 
the same momentum in neutron matter at p = 0.08 fm~"^. The temperature in all these calculations was fixed at T=5 MeV. 



CDBonn. The results for the lower energies are closer to each other. The differences in the amount of NN correlations 
is also reflected in the occupation probability n{k) defined in (fT^ . For symmetric nuclear matter at saturation density 
we obtain for n{k = 0) the values 0.89 and 0.87 for CDBonn and ArVlS, respectively. The corresponding value for 
neutron matter are n{k = 0) = 0.968 and 0.963. Earlier non-self-consistent calculations with older NN interactions 
tended to yield values of 0.83 for this quantity 3] in nuclear matter, with self-consistency raising the number to 
0.85 m. 

Examples for spectral functions A{k,u!) are displayed in Fig. |3| Again we consider the case of symmetric nuclear 
matter at the empirical saturation density and use the CDBonn interaction. As examples we consider two momenta 
k = 255 MeV/c and k = 277 MeV/c, which are below and above the Fermi-momentum kp, respectively. For the 
momentum k < kp one finds the dominant peak at an energy below the Fermi energy w = and a much smaller 
maximum at oj larger than zero. For momenta k > kp the dominant quasi-particle peak is located at positive values 
of uj and a second maximum occurs at w < 0. Since this feature of two maxima in the spectral function is reminiscent 
of the two poles that are present in the BCS approximation to the Green's function, it has been discussed as the 
formation of a pseudo-gap or as a precursor phenomenon to a pairing condensate |34j . We note, however, that our 
calculations only exhibit this feature of two pronounced maxima in the spectral function, if we consider rather low 
temperatures, in particular T < T^. The examples displayed in Fig. |31 originate from an extrapolation of the self- 
energy to r = 0. This is different from results obtained with simplified interactions, as they are used e.g. in Ref. [s^ . 
The interaction employed by Bozek yields an imaginary part of the self-energy, which is different from zero in a much 
smaller energy interval than the realistic calculations considered here. Due to this difference spectral functions with 
two maxima are obtained also at temperatures above Tc for this model interaction. 



B. Pairing correlations 

As a first step towards the study of pairing correlations, we consider the usual BCS approach. This means that 
we solve the gap equation ()23|l assuming a spectrum of sp energies e(fc) = Xfe + /x, which we determine from the 
quasi-particle energies 

e(fc) = — -f ReE(fc,e(fc) -/i) . (25) 
2m 
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FIG. 3: Spectral function for nucleons in symmetric nuclear matter at the empirical saturation density. Results of SCGF 
calculations have been extrapolated to T = 0. The CDBonn interaction has been used. 



In this equation ReE denotes the result of a SCGF calculation extrapolated to T = 0. Such spectra of quasi-particle 
energies are rather similar to the sp spectra used in other work. The corresponding BCS calculations therefore involve 
the usual procedure as it has been apphcd e.g. in Refs. @;EiEiri4, 15, 16] . 
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FIG. 4: Results for the gap-functions |Ai(fc)| for symmetric nuclear matter (p — 0.16 fm~ , upper panel) and pure neutron 
matter (p = 0.08 fm~"^, lower panel) obtained from a solution of the BCS equation H2,'i|l using the CDBonn and ArV18 
interactions at T = 0. The dotted line identifies the Fermi-momentum kp. 



Results for the gap- functions |A;(fc)| are displayed in Fig. 21 The upper panel of this figure shows results for 
symmetric nuclear matter at saturation density. The partial wave that yields the largest value for A and is therefore 
the relevant one in this case, is the '^Si — Di channel describing the proton-neutron interaction. We therefore 
display the absolute values of the gap- functions for I — and / = 2 (Aq and A2) as well as the total gap- function 
A = -^/A^ + AI as a function of the momentum k. Below we will mainly consider the value of the gap- function A at 
the Fermi momentum kp. We also compare in this figure the results obtained from CDBonn with those from ArV18. 
For smaller values of k the CDBonn yields larger values for the gap function, while ArV18 leads to larger gap values 
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for momenta larger than k = 400 MeV/c. This is true for the / = component as well as the I = 2 component and 
consequently also for the total result. This feature at large values of k is in line with our observation made above, that 
ArV18 tends to produce a larger amount of correlations at high momenta and large energies. For lower momenta, 
however, CDBonn yields larger gap- functions. Therefore the gap (at the Fermi momentum) resulting from a BCS 
calculation, which uses CDBonn (8.6 MeV) is larger than the corresponding value calculated for ArVlS (7.6 MeV), 
although ArVlS tends to produce more short-range correlations than CDBonn. 

The situation is quite similar for the case of neutron-neutron pairing in pure neutron matter, which is displayed in 
the lower part of Fig. ^ In this case the pairing effects are dominated by the ^Sq partial wave. At high momenta 
we obtain larger values for the gap function using ArVlB, whereas CDBonn yields larger values for A(fc) at low 
momenta. Therefore the value A{kp) is larger for CDBonn (1.4 MeV) than for ArV18 (1.1 MeV). These values for 
the neutron-neutron pairing gap are, however, much lower than the corresponding values for proton-neutron pairing 
at the same Fermi momentum. On one hand this sounds natural, as we know that the proton-neutron interaction is 
stronger than the neutron-neutron interaction, leading to a bound deuteron and to more correlations (see above) . On 
the other hand, however, one observes the effects of proton-proton and neutron-neutron pairing in finite nuclei, while 
there is hardly any trace of proton-neutron pairing effects in nuclei. 
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FIG. 5: The quantity Xk defined in Eq. 12611 . which represents the energy denominator for the propagator of two nucleons 
in the medium. Results are displayed for the quasi-particle approximation in the limit T — (QP, T=0), the quasi-particle 
approximation for finite temperature T = 5 MeV (QP, T=5 MeV) and the dressed propagator resulting from SCGF calculations 
(full prop., T=5 MeV). All results displayed in this figure were obtained for symmetric nuclear matter at density p = 0.16 
fm^'^, using CDBonn interaction. 



As a next step, we now try to consider the effects of temperature and short-range correlations in the solution of the 
gap equation. For that purpose we will reconsider the two-particle propagator of Eq. (|24|l but replace the spectral 
function of the superfluid phase, As(k,uj'), by the corresponding one for the normal phase, A{k,uj'). If we consider 
this propagator in the limit of the mean-field approximation (A^kjUj) = 5{lu — Sk)) at T = 0, it reduces to an energy 
denominator of the form 

1 ^"°°^A(fc,.)A(fc,.o^-^(-)-^(-') -^^-^ (26) 



-2Xfe J-oo 27r 2tt -uj - uj' -2\xk\ 

This means that the energy Xk has been defined in this equation to exhibit the effects of finite temperature and 
correlations on the two-particle propagator. Figure[Sldisplays results for this quantity Xk, the inverse of this propagator 
multiplied by -2, and compares it with \xk\ using the corresponding quasi-particle energies. The dashed-dotted line 
represents the effects of temperature, i.e. the propagator has been calculated using the quasi-particle approach for 
the spectral function and the Fermi function for the temperature under consideration, while the solid line accounts for 
finite temperature and correlation effects. The finite temperature yields an enhancement of the effective sp energy Xk 
for momenta around the Fermi-momentum, only. Including in addition the effects of correlations in the propagator. 
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we obtain larger values for Xk for all momenta. This corresponds to the well known feature that a finite temperature 
yields a depletion of the occupation probability of sp states only for momenta just below the Fermi momentum, while 
strong short-range correlation provide such a depletion for all momenta of the Fermi sea. 




% 100 200 300 400 500 

Momentum k [MeV/c] 

FIG. 6: EflFective renormalization constants for the sp strength located in the quasi-particle peak as derived from Eq. I|27|l are 
compared to the estimates of Eq. 12811 and Eq. 1291 . Further details as in Fig. |3 



Now one may try to describe this depopulation of the sp strength due to the short-range correlations by a factor 
Z, which represents the part of the sp strength contained in the quasi-particle peak. Considering the propagator of 
Eq. ^ we define 

1 , . Z'effik) 



(T, dressed) = ^^(T, QP) , (27) 

which means that we determine Z^ff e.g. from the ratio of the results displayed by the solid and the dashed-dotted 
line in Fig. El Results for such a strength factor are displayed in Fig. Typical values for Zeff{k) are around 0.8 to 
0.9 and show only a weak dependence on k. This demonstrates that the effects of correlations may be expressed in 
terms of a renormalization factor Z(k) to be used in the usual BCS equation. Such renormalization effects have been 
discussed by Bozek H^l^ and Baldo et al. 39]. Note, however, that the renormalization factor Zeff{k) is defined 
by Eq. (|27|l and the complete distribution of sp strength is required to calculate it. It will be interesting to examine, 
to which extend this factor can be approximated by simpler estimates for the strength located in the quasi-particle 
peak. 

For that purpose we present in Fig.Elalso results from simple estimates of this strength distribution. If one assumes 
that the imaginary part of the self-energy is a constant, not depending on the energy variable ui, one can estimate 
this renormalization factor by 

Zm = Y^^- (28) 

du) 

Note, however, that the real-part of the self-energy ReS calculated in SCGF yields negative slopes as a function of 
energy for various momenta and energies (see Fig.EJ. Therefore the values for Zi{k) yield values larger than 1 over a 
wide range of momenta, which is very different from the corresponding values for Zeff{k). One may try to improve 
the estimate for the strength factor by accounting for an energy dependence of the imaginary part of the self-energy 
ImE. This leads to |E2 

Z^ik) = ^ . (29) 

du; J V du J 

This improvement also does not lead to results that are consistent with Zeff{k) (see Fig. 01. This means that 
we expect results for the generalized gap-equation (|20|l which are similar to those obtained in Refs. [ssl l39l IsH 
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using appropriate values for Zi,ff{k) but we cannot give a simple reliable scheme to estimate the values for the 
renormalization constants. We therefore determine the dressed two-particle propagator in Eq. (|2(j|l within the SCGF 
approximation by extrapolating the self-energy to temperatures below Tc, as discussed above. We define an effective 
sp spectrum according to Eq. (|26|l and solve the gap equation 123f) with this effective sp spectrum. 



"1 ' I 

-- p = 0.16fm"' 

-- p = 0.08fm' 

-- p = 0.04fm"' 

p = 0.08, dressed - 

p = 0.04, dressed 
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FIG. 7: Gap parameter A{kF) in symmetric nuclear matter as a function of temperature T. Results are presented for various 
densities, with and without taking into account the dressing of the sp propagator due to short-range correlations. The pairing 
gap disappears at p=0.16 fm~^, if dressed propagators are considered. All results displayed in this figure were obtained using 
the GDBonn interaction. 



Results for the gap parameter A{kp) in symmetric nuclear matter of various densities are presented in Fig.[7|as a 
function of the temperature T. We will first discuss the results obtained within the usual BCS approximation (see 
discussion above) in the ^Si Di partial wave. At the empirical saturation density the CDBonn interaction yields 
a gap parameter A(kp) at temperature T = of 8.6 MeV (see above), which decreases with increasing temperature 
until it vanishes at T = 5.2 MeV. At p = 0.08 fm~'^, which is about half the empirical density, the value of the 
gap parameter at T = is even larger (A(/cf) = 10.6 MeV) and the gap calculated within the usual BCS approach 
disappears only at a temperature of 5.9 MeV. This increase of the pairing gap with decreasing density can be related to 
the momentum-dependence of the pairing gap A(/c) as displayed in Fig.^ The gap function increases with decreasing 
momentum. Therefore, as the Fermi momentum decreases with density, the value A(kp) tends to decrease with 
density. At even lower densities, however, this effect is more than compensated by the feature, that the phase-space of 
two-hole configurations decreases with density, so that ultimately the gap parameter will approach the binding energy 
of the deuteron in the limit of p ^ 0. This explains the decrease of the gap parameter going from p = 0.08 fm"'^ to 
p = 0.04 fm-3. 

If we take the effects of short-range correlations into account, the generalized gap equation of Eq. H20|l does not 
give a non-trivial solution for symmetrical nuclear matter at p = 0.16 fm~'^. This means that a proper treatment of 
correlation effects in nuclear matter at normal density yields a disappearance of the proton-neutron pairing predicted 
by the usual BCS approach. The effects of short-range correlations tend to decrease with density. As a consequence 
we obtain non- vanishing gaps for proton-neutron pairing at lower densities (see solid lines in Fig. 0) • This also leads 
to an increase of the critical temperature and the value of A{kp) at T = going from p = 0.08 fm^'^ to p = 0.04 fm~'^. 
Note, that these functions are qualitatively different from the corresponding BCS predictions. Differences associated 
with the various interactions are displayed in Fig. |S| For nuclear matter with a density of p = 0.08 fm~'^ the gap 
parameter is presented as a function of temperature T using the BCS approximation and the generalized gap equation 
with dressed propagators. As has already been discussed above, the ArVlS interaction yields smaller values for the 
gap parameter and the critical temperature than the CDBonn interaction. 

Effects of pairing correlations on the spectral function are visualized in Fig. |51 As an example we consider nuclear 
matter at p = 0.08 fm~'^ and show results for the spectral function without {A(k,u!)) and with inclusion of pairing 
correlations {As{k,u!), see Eq. The momentum considered for this figure, k = 193 MeV/c, is slightly below the 

Fermi momentum kp = 208 MeV/c. One observes that the inclusion of pairing correlations enhances the maximum 
of the spectral distribution at positive values of lu considerably and shifts the quasi-particle peak to more negative 
values of lu. The pairing correlations modify the spectral distribution into the direction which is obtained in the simple 
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FIG. 8: Gap parameter A(fcF) in symmetric nuclear matter at p=0.08 fm^'^ as a function of temperature T using the ArV18 
and the CDBonn interaction. 
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FIG. 9: Spectral function for nucleons with momentum k — 193 MeV/c with (solid line) and without (dashed line) inclusion 
of pairing correlations. Results are presented for nuclear matter of p = 0.08 fm~"^ at a temperature T = 0.5 MeV. 



BCS approximation for As{k,uj) in Eq. (|21|l . Also note, that these modifications of the spectral function As{k,uj) 
as compared to A{k,Lu) is limited to a small interval of energies around lu — and to momenta close to the Fermi 
momentum. 

Finally, we consider the case of neutron-neutron pairing in pure neutron matter. We will focus the attention to 
densities, where the pairing correlations in the ^5*0 partial wave are dominating. Results for the gap parameter Aikp) 
as a function of temperature are displayed in Fig. llUI Using the BCS approximation with sp energies derived from the 
quasi-particle energies of SCGF calculations we obtain a gap at T = 0, which, for the range of densities considered, 
increases with decreasing density. This is in agreement with results of similar calculations, which are summarized e.g. 
in 0. 

The effects of short-range correlations are weaker in neutron matter than in nuclear matter. This has been discussed 
already above in connection with the results displayed in Fig. [3 This can also be seen in a comparison of the dressed 
two-particle propagator and the effective strength factor Zeff defined in Eq. (|27|l . While a calculation of Z^/f for 
symmetric nuclear matter at p = 0.16 fm^"^ yields values for Zf,ff, which are typically around 0.8 (see Fig. 1^, 
corresponding values for Z^ff in neutron matter at p = 0.08 fm~*^ are around 0.9. Nevertheless, also these weaker 
effects of short-range correlations in neutron matter are sufficient to suppress the formation of a pairing gap in neutron 
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FIG. 10: Gap parameter A(A;f) in neutron matter as a function of temperature T. Results are presented for the usual BCS 
approximation (dashed lines) and the solution of the generalized gap equation i'ZUl in the ^So partial wave using the CDBonn 
interaction. 



matter at p = 0.08 fm~^. Such a suppression of pairing correlations is also observed at smaller densities. In this case, 
however, the inclusion of the correlation effects just leads to a reduction of the gap parameter at a given temperature 
and a reduction of the critical temperature (see Fig. I10|l . Extrapolating our results to neutron matter with higher 
densities, we expect that the short-range correlations will suppress the formation of pairing in the F2 partial 

waves at those densities [s^. 



V. CONCLUSIONS 



An attempt has been made to treat the effects of short-range and pairing correlations in a consistent way within the 
T-matrix approach of the self-consistent Green's function (SCGF) method. The pairing effects are determined from 
a generalized gap equation that employs sp propagators fully dressed by short-range and tensor correlations. This 
equation is directly linked to the homogeneous solution of the T-matrix equation of NN scattering in the medium, 
which is one of the basic equations of the SCGF approach at temperatures above the critical temperature for a phase 
transition to pairing condensation. While short-range and tensor correlations yield a redistribution of sp strength 
over a wide range of energies, the effects of pairing correlations on the spectral function are limited in nuclear matter 
to a relatively small interval in energy and momentum around the Fermi surface. 

The formation of a pairing gap is very sensitive to the quasi-particle energies and strength distribution at the Fermi 
surface and can be suppressed by moderate temperatures. The formation of short-range correlations are sensitive to 
a larger range of energies and momenta. So we observe, that the non-local CDBonn interaction is softer with respect 
to the formation of short-range correlations but yields larger pairing gaps compared to the local ArV18 model for the 
NN interaction. 

From this sensitivity to different areas in momentum and energy one may conclude that the features of short- 
range correlations should be rather similar in studies of nuclear matter and finite nuclei. The investigation of pairing 
phenomena, however, is rather sensitive e.g. to the energy spectrum around the Fermi energy. Therefore the shell 
effects of finite nuclei may lead to quite different results for pairing properties than corresponding studies in infinite 
matter. 

The redistribution of sp strength due to the short-range correlations has a significant effect on the formation of a 
pairing gap. While the usual BCS approach predicts a gap for proton-neutron pairing in nuclear matter at saturation 
density as large as 8 MeV, the inclusion of short-range correlations suppresses this gap completely. Correlation effects 
are weaker at smaller densities, but still lead to a significant quenching of the proton-neutron pairing gap and to a 
reduction of the critical temperature for the phase transition. Compared to symmetric nuclear matter correlation 
effects are weaker in neutron matter. Nevertheless, the inclusion of correlations suppresses the formation of a gap for 
neutron-neutron pairing at p = 0.08 fm^'^ completely and yields a significant quenching at lower densities. 
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The effects of dressed sp propagator in the generahzed gap equation could be described in terms of an effective 
strength factor Z^ffk, which has been considered in the hterature before |3^|3^|^. Unfortunately, we have not 
been able to derive the value of this strength factor from bulk properties of the self-energy. 

Although the effects of pairing correlations on the sp Green's function is weak and limited to a small range in energy 
and momentum, these modifications are very important to extend SCGF calculations to densities and temperatures 
that suffer from the so-called pairing instability. The present study is a first step towards a consistent treatment of 
pairing and short-range correlations. 

This work is supported by the U.S. National Science Foundation under Grant No. PHY-0140316 and the "Landes- 
forschungsschwerpunkt Quasiteilchen" of the state of Baden Wiirttemberg. H.M. would like for the hospitality of the 
Department of Physics of Washington University in St. Louis, where a major part of this work has been done. 
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